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ABSTRACT 

The number of low-mass brown dwarfs and even free floating planetary mass objects in young 
nearby star-forming regions and associations is continuously increasing, offering the possibility to 
study the low-mass end of the IMF in greater detail. In this paper, we present six new candidates 
for (very) low-mass objects in the Taurus star-forming region one of which was recently discovered 
in parallel by Luhman et al. (2009b). The underlying data we use is part of a new database from a 
deep near-infrared survey at the Calar Alto observatory. The survey is more than four magnitudes 
deeper than the 2MASS survey and covers currently ~1.5 deg^. Complementary optical photometry 
from SDSS were available for roughly 1.0 deg^. After selection of the candidates using different color 
indices, additional photometry from Spitzer/IRAC was included in the analysis. In greater detail we 
focus on two very faint objects for which we obtained J-band spectra. Based on comparison with 
reference spectra we derive a spectral type of L2±0.5 for one object, making it the object with the 
latest spectral type in Taurus known today. From models we find the effective temperature to be 
2080±140 K and the mass 5-15 Jupiter masses. For the second source the J-band spectrum does not 
provide a definite proof of the young, low-mass nature of the object as the expected steep water vapor 
absorption at 1.33 /im is not present in the data. We discuss the probability that this object might 
be a background giant or carbon star. If it were a young Taurus member, however, a comparison to 
theoretical models suggests that it lies close to or even below the deuterium burning limit (<13 Mj„p) 
as well. A first proper motion analysis for both objects shows that they are good candidates for being 
Taurus members. 

Subject headings: stars: low-mass, brown dwarfs, stars: pre-main sequence, stars: formation, planetary 
systems: formation, planetary systems: protoplanetary disks 



1. INTRODUCTION 

Taurus is certainly one of the best studied low-mass 
star-forming (SF) regions in the northern hemisphere. It 
lies at a distance of ~ 140 pc (Kenyon et al. 1994) and 
is very young (median age ~1 Myr, Briceho et al. 2003; 
Luhman et al. 2003). The initial mass function (IMF) in 
Taurus has been studied intensively in the last years us- 
ing large surveys mostly in the optical. It first seemed as 
if at the low-mass end the number of brown dwarfs (BDs) 
was too small compared to the expectations for a typical 
IMF (Briceiio et al. 2002). While other studies found in- 
Electronic address: quanz@astro.phys.othz.ch 

^ Based on observations made at the Calar Alto Observatory 
Based on observations made with ESO Telescopes at the Paranal 
Observatories under program ID 278.C-5043A. This work is based 
in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. 



dications that this might not be the case (Luhman et al. 
2006a; Guieu et al. 2006) a recent analysis finds that a 
surplus of late K and early M type stars seems still to 
be present (Luhman et al. 2009a). Apart from one ex- 
ception the least massive objects in Taurus known today 
have spectral types of M9-M9.5 with corresponding ef- 
fective temperatures of ~ 2400-2500 K (Luhman et al. 
2009b; Luhman 2006; Guieu et al. 2006, and references 
therein). Table Al shows that theoretical models predict 
these objects to have masses above the deuterium burn- 
ing limit of ~0.013 Mo«13 Mj„p (Burrows et al. 1997) 
which makes them very low-mass brown dwarfs. Only 
very recently Luhman et al. (2009a) reported the detec- 
tion of the first LO dwarf in Taurus where the spectral 
type was derived from an I-band spectrum. However, 
this object appears to be underluminous for its assumed 
young age and we will discuss this object and its proper- 
ties in more detail in section 4. In theory, young objects 



2 



Quanz et al. 



(1-3 Myr) with L and T spectral types and corresponding 
effective temperatures below ^2300 K could potentially 
be of planetary mass (Table A2). 

For other young star-forming regions the detection of 
L- and T-type objects has already been reported and 
their numbers are rather small but continuously increas- 
ing. In the 3-5 Myr old cluster cr-Orionis one object, S 
Ori 70, has an observed spectral type of T5.5±l (Zap- 
atero Osorio et al. 2002), one object is of spectral type 
L4 and three objects have spectroscopically confirmed 
spectral types of L0-L1.5 (Zapatero Osorio et al. 2000). 
If these objects are members of the young cluster then 
they can be considered as very good candidates for iso- 
lated planetary mass objects (i.e., with masses M<13 
Mjup). However, McGovern et al. (2004) found for one 
potential candidate in cr-Orionis (S Ori 47) spectroscopic 
evidence that the object is rather an object from the 
field than from the young cluster. Recently, Bihain et al. 
(2009) report the detection of three additional objects 
in cr-Orionis (S Ori 72-74) with theoretical masses of a 
few Jupiter masses (with L-T spectral types), but spec- 
troscopic confirmation is still pending. In addition to 
the objects in cr-Orionis Lodieu et al. (2008) and Luh- 
man et al. (2008) report the finding of L-type objects in 
the Upper Sco Association and the Chamck;on I region, 
respectively (median ages ~3-5 Myr). Weights et al. 
(2009) report the spectroscopic confirmation of 38 very 
low mass objects below the hydrogen burning limit in the 
Orion Nebular cluster (median age ^1 Myr). 10 of their 
objects could have masses below the deuterium burn- 
ing limit and be thus planetary mass objects. Burgess 
et al. (2009) recently published the detection of a T-dwarf 
candidate in the 3 Myr-old star-forming region IC 348 
based on methane-band imagery. The found a prelim- 
inary spectral type of T6 for the object and concluded 
that the detection of this source is consistent with the 
extrapolation of current lognormal IMF estimates down 
to the planetary mass domain given the survey area and 
completeness of their sample. Furthermore, there are 
three optically classified L dwarf candidates in the field 
showing some spectral signatures of youth (Kirkpatrick 
et al. 2001, 2006; Luhman et al. 2006b). Four addi- 
tional objects in young star-forming regions (1-3 Myr) 
are at the borderline between M and L-type objects but 
are subject to controversial discussions in the literature: 
Jayawardhana & Ivanov (2006a, b) derived for four ob- 
jects in Chameleon, Lupus and Ophiuchus spectral types 
of LO and planetary masses while Luhman et al. (2007) 
and AUers et al. (2007) derived at least for three of the 
objects significantly earlier spectral types (M7.25-M8.75) 
and higher masses. In the and the TW Hya association 
(5-10 Myr) only late M-typc and thus more massive ob- 
jects were found (Liu et al. 2003; Gizis 2002). 

In this paper we present new candidates for low-mass 
BD in Taurus. At least two objects are good candidates 
for being among the least massive objects known today 
in this SF region possibly lying below or close to the 
edge of isolated planetary mass objects. For one object 
we derive a spectral type of L2±0.5 and a most likely 
mass between 5 and 15 Jupiter masses. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. OMEGA2000 data 



In a new deep near-infrared (NIR) survey of Taurus 
we focused on regions of comparably high extinction 
where previous optical surveys might have missed po- 
tentially young, very low-mass objects. In the first part 
of our ongoing survey we observed the vicinity of 55 
dense molecular cores (Onishi et al. 2002) between Oc- 
tober 2004 and January 2006 with the wide-field NIR 
camera Omega2000 on the 3.5 meter telescope in Calar 
Alto (Spain)^ . The camera is equipped with a 2048 x 2048 
pixel HAWAII-2 detector, has a pixel scale of 0.45"/pixel 
and an effective field-of-view of ~15'. The 55 dense cores 
were covered with 24 telescope pointings. One additional 
field harboring a proto-brown dwarf (Apai et al. 2005) 
was included in the survey. Up to now the NIR observa- 
tions cover an area of ~1.5 deg^ on the sky (Table A3). 
Each field was observed in the three NIR filters J, H, and 
Kg with central wavelengths at 1.21, 1.65, and 2.15/im, 
respectively. For two fields the Kg-hand data had to be 
disregarded in the data analysis as cloudy skies disturbed 
the observations significantly. Following a pre-defined 
dither pattern 30 frames were taken on each field and 
filter. Each frame had an on-source integration time of 
one minute, so that after co-addition the total integra- 
tion time amounted to 30 minutes per field and filter. 
The limiting magnitudes were ~22.0 mag in J, ~20.5 
mag in H and ~ 19.5 mag in Kg. The depth of the sur- 
vey allows us to detect young objects (1-5 Myr) down 
to a few Jupiter masses in all three filters if they were 
present. 

The data reduction was carried out with a dedicated 
pipeline developed at the MPIA. The pipeline runs 
within the MIDAS environment and does flat-fielding, bad 
pixel and cosmic ray removal, sky subtraction and finally 
weighting and co-addition of the individual images to a 
combined final image for each field and filter. If one 
frame suffers substantially from bad quality, e.g., due 
to bad observing conditions in terms of seeing or trans- 
parency, this frame is rejected and not included in the 
final image. A more detailed overview concerning the 
pipeline is provided in the Omega2000 Handbook^. 

Astrometry was applied to each field by comparing the 
positions of detected 2MASS point sources to those listed 
in the 2MASS Point Source Catalog (PSC) (Cutri et al. 
2003). On average the precision of the astrometry is 
better than 0.25" with respect to the known 2MASS po- 
sitions. 

After visual inspection of the resulting final frames, 
photometry was carried out with the daophot package 
within the IRAF^ environment. As the point spread func- 
tion (PSF) does not vary significantly over the chip, PSF- 
photometry was carried out instead of simple aperture 
photometry, in order to ensure reasonable fiux measure- 
ments of close binaries and in more crowded regions. A 
reference PSF was created for each filter and field indi- 
vidually and then applied to all sources detected with a 
3-cr confidence level. The resulting tables containing the 
fluxes for each fleld and fllter were then cross-matched in 
order to identify sources that were detected in all three 
filters. For the cross-matching the coordinates of the 
measured flux peaks were required to match within 2 

2 http://w3.caha.es/CAHA/Instruments/O2000/index2.html 
^ http://www.caha.es/CAHA/Instruments/O2000/index.html 
^ http;//iraf.noao.edu/ 
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pixels. 

Finally, the photometric calibration was achieved by 
comparing the measured instrumental magnitudes to 
those listed in the 2MASS PSC. Here, only objects with 
the best "quality flag" and "read flag" in the 2MASS PSC 
were considered serving as photometric reference sources. 
Typically between 50 and 100 objects satisfying these cri- 
teria were detected in each Taurus field. By plotting the 
difference between the observed instrumental magnitudes 
and those measured in 2MASS in one filter against an 
observed color (e.g., Jo2K - J2MASS vs. J02K - Ho2k) 
the color correction terms between the Omega2000 and 
2MASS filter systems were derived. This procedure was 
done for each observed Taurus field individually as ob- 
serving conditions did vary between the individual fields. 
The mean transformation errors for each field and filter 
are also given in Table A3. 

2.2. SDSS data 

To differentiate between potential young members 
of Taurus and background and foreground contami- 
nants/sources it is crucial to rely not only on NIR 
data but to take into account also optical informa- 
tion. The Sloan Digital Sky Survey (SDSS) offers up to 
five additional bands in the optical wavelength regime^ 
(Finkbeiner et al. 2004) and covers parts of our NIR sur- 
vey area. In the upper part of Table A3 we provide 
estimates for the overlap between SDSS and our fields 
(in total ~1.0 deg^). For the rest of the paper we will 
only focus on the overlapping regions. We matched our 
Omega2000 NIR catalog with the SDSS data allowing 
a maximum offset of 2" for an objects' position in both 
data sets. 

2.3. Spitzer IRAC data 

Photometry at longer IR wavelengths allows us to de- 
tect the presence of disks around potential Taurus can- 
didates supporting a young age of the object and thus 
membership of the SF region. Images taken with the 
IRAC instrument onboard the Spitzer Space Tele- 
scope were downloaded from the data archive^' for the 
finally selected objects (see section 3). The so-called 
"Post Basic Calibrated Data (PBCD)" were used to de- 
termine the flux of the individual candidates in the first 
two IRAC filters at 3.6 and 4.5 /xm. Unfortunately, the 
sensitivity of these images is insufficient to detect some 
of the fainter candidates at 5.8 and 8.0 /im. The initial 
pixel values of the images were converted from MJysr"^ 
to count rates DNs~^ (Data Numbers per second). Af- 
terwards, aperture photometry was carried out using the 
standard atv.pro routine in IDL with an aperture size 
of 2 pixel and a sky annulus from 10-20 pixel. Aper- 
ture corrections were carried out as described in the 
IRAC Data Handbook to obtain the final magnitudes via 
m = —2.5 log(a;) -|- Azp with x denoting the fiux mea- 
sured in DNs~^ and Azp being the zero point for each 
filter^. 

^ http:/ /photo. astro. princcton.edu/orioiidatarclcasc/ 

® These data arc part of a large scale IRAC Taurus survey (PI: 

Deborah Padgett) and publicly available. 

Zero points taken from Hartmann et al. (2005); 19.66 (3.6 /im), 

18.94 (4.5/im), 16.88 (5.8/im), 17.39 (8//m). 



2.4. VLT/ISAAC data 

For two of the three least luminous, and hence most 

interesting, objects we were granted Director's Discre- 
tionary Time (DDT) for ISAAC on ESO's VLT to obtain 
J-band low-resolution spectra. In low-resolution mode, 
ISAAC has a pixel scale of 0.147 "/pixel and a slit length 
120". The observations took place between February 7 
and February 11, 2007. The spectra cover a wavelength 
range between 1.1 and 1.4 /zm and we chose a slit width 
of 0.8" resulting in a resolution of i? « 680. The observa- 
tions were done in a nodding mode, with a 40" nodding- 
throw along the slit. The detector integration time (DIT) 
was 120 s for both objects. The number of DIT expo- 
sures averaged out within the same frame (NDIT) and 
the number of frames (NINT) were 2 and 10, and 2 and 6 
for CAHA Tau 1 and CAHA Tau 2, respectively The air- 
mass was between 1.5 and 2.0 in general, but increased 
up to 2.3 for the observations of CAHA Tau 2. The 
seeing was mostly better than 0.8". The data reduc- 
tion was done with the IRAF package one spec following 
the standard procedure for low-resolution spectra: sky 
subtraction, fiatfielding, cosmic ray removal, and finally 
extraction of the spectrum using an interactively fitted 
aperture. After wavelength calibration using argon lines 
the individual spectra for each object were averaged in 
combination with a 3-sigma clipping for the individual 
points. Telluric lines were corrected afterwards using 
standard stars^ that were observed roughly at the same 
airmass directly after the science objects. To correct for 
the spectral response of the detector we compared the 
observed spectra of the standard stars with theoretical 
models of the same spectral type and applied correction 
factors to obtain an identical spectral slope. 

3. RESULTS 

3.1. Candidate selection 

To identify potential low-mass objects of the Taurus 
region among the thousands of detected objects an op- 
tical/NIR color-magnitude diagram (CMD) served as 
starting point. In particular for fainter objects one has 
to find a way to distinguish between nearby low-mass ob- 
jects, distant giant stars and also distant galaxies. While 
even distant galaxies with high red-shifts typically have 
I — J colors between and 2.5, nearby low-mass ob- 
jects normally appear significantly redder. Figure Al, 
where we plot K against I — J, takes advantage of this 
fact. As the overplotted theoretical isochrones (see be- 
low) are provided in the CIT photometric system, the 
NIR data in this figure were transformed into the CIT 
system. The optical data is in the Johnson-Cousin sys- 
tem. The required equations to transform between the 
different photometric systems are given in the appendix. 
In the plot we show only objects with a photometric error 
< 0.1 mag in the NIR and with a S/N > 5 in the SDSS 
i and z band. Out of more than 9000 objects that were 
commonly detected in the SDSS database and the NIR 
observations 5251 sources fulfilled these criteria. Theo- 
retical isochrones computed by Baraffe et al. (2002) and 
Chabrier et al. (2000) are overplotted in the CMD. The 
isochrones have an age of 5 Myr, are scaled to a distance 
of 140 pc for Taurus (Kenyon et al. 1994) and cover the 

* We used the B5 stars Hip023151 and Hip023946. 
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mass regimes between 0.02 Mq and 1.4 Mq (BarafFe et al. 

2002) , and 0.002 M© and O.OTSAf© (Chabrier et al. 2000). 
The assumed age is a rather conservative assumption be- 
cause, as mentioned in the introduction, Taurus seems 
to be significantly younger (r ~ 1 Myr, Briceno et al. 

2003) . Since younger objects tend to appear brighter 
(they have larger radii) and redder (e.g., due to circum- 
stellar material), objects falling along the isochroncs or 
lying above them are potential Taurus candidates and 
should be selected. As expected the least luminous ob- 
jects are not selected as they appear too blue in / — J 
and arc thus background stars or galaxies. In total, using 
the isochroncs as a first selection criterion, the number 
of objects was reduced from 5251 down to 646. 

These objects are now, in a second step, plotted in a 
NIR color-color diagram (Figure A2). Also here the CIT 
photometric system is used to be consistent with the first 
selection step. It can be seen that most of the selected 
objects lie along the reddened main sequence or along 
the reddened giant branch. Only very few sources pop- 
ulate the classical TTauri star locus (Meyer et al. 1997) 
or show otherwise any NIR excess emission. This indi- 
cates that the sample seems to be mostly contaminated 
by (background) stars. To minimize the possibility of se- 
lecting these objects, we considered only objects lying on 
the right hand side of a reddened M3 main sequence star 
and showing positive NIR colors. This selection resulted 
in 22 objects which were then subject to an individual 
inspection on the images. As some objects were too close 
to the edge of the detector in some filters and one object 
was observed twice in two adjacent telescope pointings 
the final sample of candidates consisted of 16 objects. 
Their positions in the color-color-diagram are indicated 
by crosses and squares in Figure A2. 

Finally, the coordinates for the 16 remaining ob- 
jects were crosschecked with the SIMBAD Astronomi- 
cal Database^. 10 objects were already listed (boxes in 
Figure A2), but 6 candidates did initially not have a 
counterpart in the database (crosses in Figure A2). 

For the 10 previously known objects the names, mag- 
nitudes, spectral types and references are summarized in 
Table A4 along with our new NIR photometry. We stress 
that the objects with known spectral type are without ex- 
ception low-mass stars and BDs in Taurus. This empha- 
sizes the quality of our selection criteria. With KNPO- 
Tau 4 we even re-detected an object with a spectral type 
as late as M9.5. However, since we applied rather strict 
selection criteria, we did not re-identify all of the known 
low-mass objects and BDs that lie in our observing fields. 
Some of them have for instance bluer H — K colors than 
allowed by our selection criteria. One example can be 
seen in Figure A3. In Table A4 also the magnitudes 
from the 2MASS Point Source Catalogue are given for 
comparison. One can see that some objects appear quite 
variable while other sources apparently did not change in 
brightness at all. One possible explanation for the vari- 
ability of the objects could be strong accretion events. 
CFHT-BD-Tau 19, for instance, shows clear evidence for 
accretion in its optical spectrum (Guieu et al. 2006). 

For the 6 objects that were initially not listed in SIM- 
BAD the coordinates and photometric properties are 
shown in Table A5 together with estimates for their 

^ http: / / simbad.u-strasbg.fr / simbad / 



optical extinction based on extinction maps. However, 
during the revision process of this paper Luhman et al. 
(2009b) also published results for our candidate CAHA 
Tau 6 and found the object to be very young BD with a 
spectral type of M9.25 and a mass of ~0.015 M©. Most 
interestingly, the object turned out to be a member of an 
isolated BD binary system with FU Tau (now FU Tau 
A) being the primary object which we also re-identified 
among the 10 previously known objects (see Table A4). 
Although CAHA Tau 6 (now FU Tau B) was thus al- 
ready confirmed to be a young low-mass object, we keep 
it in our sample in the following, especially because no 
NIR data has been published so far. Thus, the properties 
of all our initially selected six candidates are discussed 
in greater detail in the following sections. NIR finding 
charts for CAHA Tau 1-5 are available from the online 
material. In Luhman et al. (2009b) an optical image of 
the CAHA Tau 6 region is provided. 

3.2. Photometric properties of candidates 

In the following we will discuss the photometric prop- 
erties our selected candidates. From now on we use the 
2MASS photometric system in the NIR throughout the 
rest of the paper as most reference sources were pub- 
lished with 2MASS NIR data. In Figure A3 we show 
again a NIR color-color diagram. This time, however, we 
compare the colors of our six candidates to those of pre- 
viously detected BDs in Taurus from Guieu et al. (2006) 
and embedded carbon stars from Liebert et al. (2000). 
As we are interested in very low-mass objects we show 
only objects from Guieu et al. (2006) with spectral types 
of M7 or later. We notice two things: First, the colors 
of most of our candidates do agree quite well with low- 
mass objects. Only CAHA Tau 5 appears a little too 
blue in J — H. Second, some of the carbon stars have 
NIR colors that are very similar to those of our young 
low-mass candidates. Especially CAHA Tau 3 (and to a 
lesser extend also CAHA Tau 4) finds itself surrounded 
by carbon stars. 

In Figure A4 we continue the comparison of our can- 
didates to other young low-mass sources. This time, we 
increase the wavelength range in the NIR by using IRAC 
data at 3.6 and 4.5 /^m and we include optical colors in 
the diagrams. Furthermore, we add young, planetary 
mass candidates in a Orionis from Zapatero Osorio et al. 
(2007) as reference objects. It shows that in the left-hand 
plot of Figure A4 our candidates populate nicely the 
same region of the color-color diagram as the other young 
low-mass objects. Only CAHA Tau 6 appears much red- 
der in the IRAC color but Luhman et al. (2009b) argued 
that this object is indeed particularly young (~ 1 Myr) 
and shows significant infrared emission excess. In the 
right-hand plot, where we have chosen a broader wave- 
length baseline for the colors, our candidates tend to oc- 
cupy the left region of the plot, i.e., they are relatively 
blue in Kg ~ 4.5 /im compared to the reference objects. 
However, the colors are in general agreement with those 
of the other young low-mass objects. Only CAHA Tau 4 
lies a little offsite compared to the other candidates and, 
again, CAHA Tau 6 appears very red in the Ks — 4.5 /xm 
color. 

3.3. J -hand spectroscopy 
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As seen for instance in AUcrs et al. (2006, 2007) one 
has to be careful if one tries to derive stellar param- 
eters of young objects such as spectral types, effective 
temperatures and masses solely from photometric mea- 
surements. Hence, we complemented our data with NIR 
J-band spectra for two of the three least luminous ob- 
jects. We focused on these objects as they are potentially 
also the least massive and hence, for our purpose, most 
interesting candidates. 

In Figures A5 we show the J-band spectra of CAHA 
Tau 1 and CAHA Tau 2 and compare them to data of 
field BDs and giants (left), and to young objects pre- 
sented by AUers et al. (2007) and model spectra from 
Burrows et al. (2003, 2006) (right). The left hand side 
of Figure A5 suggests that our candidates are late-type 
objects, as the spectrum of the early M-type giant HD 
120052 does not show the deep water absorption band at 
the red edge of the spectrum that is seen in our candi- 
dates. Also, the general slope of the spectra is in agree- 
ment with late type objects. Furthermore it seems clear 
that our candidates have a low surface gravity, an indi- 
cation for a possibly young age, as no strong spectro- 
scopic features are present. However, in the spectrum of 
CAHA Tau 1 we find weak spectral signatures: There 
is a dip near the position of the Na I doublet around 
1.14 /xm and also hints for the K I doublets near 1.17 - 
1.18 /im and 1.24 - 1.25 /im which show up much stronger 
in the spectra of the higher surface gravity field stars. In 
Figure A6 a direct comparison between the J-band spec- 
trum of CAHA Tau 1 and those of late type field objects 
is shown. The overall very similar shape and slope is 
striking and also the existence of similar spectral fea- 
tures, even if they are very weak in case of CAHA Tau 1. 
The best overall agreement in shape and slope between 
a reference spectrum and CAHA Tau 1 is found for Kelu 
1 AB. Based on this we estimate the spectral type of 
CAHA Tau 1 to be L2±0.5 which is the latest spectral 
type known in Taurus today. We will discuss this finding 
in more detail in section 4.1. 

A spectral classification for CAHA Tau 2 is difiicult as 
its spectrum lacks any spectral lines. However, it is not 
unusual for young low-mass sources to have almost no 
spectroscopic features as can be seen in McGovern et al. 
(2004) and Lodieu et al. (2008). The first group pub- 
lished spectra of the young Taurus object KPNO-Tau 4 
and of (J Ori 51 while the second group presented spec- 
tra of young L dwarfs in Upper Sco. All of these objects 
showed barely any spectral features in the J-band, pre- 
sumably due to their low surface gravity, and appear in 
this respect similar to CAHA Tau 2^°. Finally, a com- 
parison to theoretical models and reference spectra of 
young low mass objects further emphasizes the almost 
featureless behavior of these objects (right hand side of 
Figure A5). 

One thing, however, appears strange in the spectrum 
of CAHA Tau 2 and so far we lack a satisfying expla- 
nation if it was a young object. The drop in the flux 
caused by water vapor around 1.33 ^m is not as steep 
as in the other low- mass objects including CAHA Tau 1. 
Although the peak fiux in the J-band is reached between 

^'^ As mentioned above we re-identified KPNO-Tau 4, a con- 
firmed young and low-mass brown dwarf, during our selection pro- 
cess (see Table A4). 



1.32 and 1.33 /im, as it is also the case for the reference 
BDs, the flux of our candidate decreases more slowly with 
increasing wavelength. This behavior is not typical for 
young, low- mass objects (see, e.g., Lodieu et al. 2008). 
However, J-band spectra of C-N and C-J carbon stars as 
presented by Tanaka et al. (2007) show a similar behavior 
with no significant drop in the flux up to 1.35 /im. And 
also the flux of the M7 HI star HD 108849 in Figure A5 
shows a rather modest decline longward of 1.33 /um in 
addition to an overall featureless J-band spectrum. 

To summarize, for CAHA Tau 1 we find spectroscopic 
support for the young and low-mass nature of the object 
from both, the weak spectral features and the overall 
slope of the spectrum, and we can estimate a spectral 
type. For CAHA Tau 2 a definite spectroscopic proof 
that the object is indeed a young, low-mass Taurus mem- 
ber is pending and although the general slope is very 
similar to that of CAHA Tau 1 the lack of strong water 
absorption around 1.33 /um talks in favor of a background 
carbon star or giant. 

3.4. Proper motion of Taurus members and our 
candidates 

Another way to find additional support for the low- 
mass nature of our candidates is to analyze their proper 
motion. For CAHA Tau 6 Luhman et al. (2009b) already 
showed that it most likely shares common proper motion 
with FU Tau A which itself shows a proper motion con- 
sistent with Taurus membership. Thus we focus on the 
remaining 5 objects. 

Proper motion studies of stellar members of the 
whole Taurus- Auriga association show a wide disper- 
sion in proper motions along the right ascension axis, 
with — 5 < < 15mas/yr, but a consistently nega- 
tive proper motion along the declination axis, with 
—35 < iJ,s < — 5mas/yr (Frink et al. 1997), mostly due 
to the reflex motion (Jones & Herbig 1979). We have 
run a simulation of our flelds using the Besangon model 
(Robin et al. 2003, 2004), adding an extinction layer at 
150 pc distance with variable (extinctions. The model 
does not include the star forming region, only the smooth 
components of the Galactic disks and halo. This simula- 
tion (see Figure AlO) shows that most giant and sub- 
giant stars have proper motion < 5mas/yr, cor- 
responding to halo-type velocities (at most) and dis- 
tances larger than lOkpc (see subsection 4.2). The frac- 
tion of giant and sub-giant stars with proper motion 
^15 < -lOmas/yr and 16 < J < 18 is 3 ± 1% (l-a 
statistical error). The fraction for the brighter M gi- 
ants would be even smaller. Therefore, a significantly 
larger proper motion measurement would reject the gi- 
ant hypothesis, while a proper motion fig > — 8mas/yr 
or lual > 30mas/yr would most likely exclude a Taurus 
membership. Because the Taurus proper motion space 
overlaps with that of disk dwarfs, proper motions cannot 
completely remove the field dwarf contamination. A frac- 
tion of 5% of field dwarfs with 16 < J < 18 have a proper 
motion of < 30mas/yr and /i^ < — 15mas/yr. 

To complement our near-infrared and SDSS images, 
we used public data from CFHT's 12k and MegaCam 
cameras, and, for CAHA Tau 1, the catalogues from 
the UKIDSS Galactic Plan Survey K band data from 
the Data Release 3 (Lawrence et al. 2007). CFHT im- 
ages were detrended using the Elixir pipeline (Magnier & 
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Cuillandre 2004). Object detection and astrometry was 
performed using SExtractor (Bertin & Arnouts 1996) on 
all images (CFHT, SDSS, Omega2000). Bright objects 
(SNR>100, a few dozens) were used to derive a prelim- 
inary transformation, which is refined with all objects 
with SNR>10 (one hundred). For the three faintest tar- 
gets, 2MASS data have too large positional uncertainties 
to improve the fit. The photonoise errors are negligible 
at our targets' brightness; the error budget is dominated 
by transformation errors, which includes refraction, in- 
strumental distortions and proper motions of the refer- 
ence stars. Therefore we set the positional accuracies to 
a minimum of 30 mas, typical of the best position dis- 
persion measured by scamp for bright stars. This dis- 
persion rises to 80 mas in some cases (e.g. when using 
SDSS as the reference catalogue). We then fit a proper 
motion on the CFHT, SDSS, O2000 and UKIDSS po- 
sitions when a detection was found within 1 arcsec of 
the O2000 position, for all objects with SDSS z and 
O2000 Kg magnitudes similar to our candidates (see Ta- 
ble A6 and Figure All). We have not corrected for the 
velocities of the stars which form our reference grid and 
we present relative proper motions only. The Besangon 
simulation shows that stars brighter than J < 16 mag 
have a small averaged proper motion of {Ha) = 2 mas/yr 
and (ns) = — 4mas/yr, although some have large proper 
motions. 

The data set used for these measurements is far from 
optimal: we compare positions obtained with four differ- 
ent instruments, with the first epochs obtained in the red 
and the last epoch(s) obtained in the near- infrared. We 
did not find any correlation between positional shifts and 
hour angles, or filter bands, so that most of the transfor- 
mation error budget is likely to be due to the imperfect 
corrections of the (relative) instrumental optical distor- 
tions, and the proper motion dispersions of the refer- 
ence objects. Because we could not calculate a global 
astrometric solution simultaneously for all instruments, 
proper motion measurements could not be determined 
during the registration process itself. 

3.4.1. CAHA Tau 1 and 2 

Despite these limitations we find that the proper mo- 
tions of CAHA Tau 1 and CAHA Tau 2 are compati- 
ble with typical Taurus members' proper motions (Fig- 
ures AlO, All). In order to judge whether the giant 
hypothesis is rejected, we need to know how significantly 
they differ from ji = 5 mas/yr. Because of the short time 
baseline and the required proper motion accuracy, it is 
important to have a good estimate on the proper motion 
errors, which will be smaller than the signal by a factor 
of a few at most. A first estimate of the positional accu- 
racy of each data set is given by the shift dispersion of 
the bright reference stars, subtracting the estimated ref- 
erence catalogue errors, to which we add the photonoise 
of the candidate. However, this might be a lower limit 
for CAHA Tau 1 as it is located in a corner of the CFHT 
CCDs, in an empty area. The location of CAHA Tau 2 
is more central and favorable. The colors of the candi- 
dates are typical of the bulk of the reference stars. Many 
reference stars are background stars with significant red- 
dening, and therefore red colors. Furthermore, the dif- 
ference in colors does not correlate with proper motions 
or dispersions around the proper motion fit. The large 



number of CFHT images at a given epoch for CAHA 
Tau 1 allows us to measure the positional dispersion. As 
described above, we used this value as a proxy for the 
positional accuracy. Unfortunately, this robust method 
cannot be used for Omega2000, and even less for the 
UKIDSS single image. Moreover, the UKIDSS DR3 cat- 
alogue comes with no astrometric errors, and the proper 
motion fit relies crucially on this single point, as it ex- 
tends by one year the time baseline. We derive the proper 
motion error from the dispersion of the proper motions 
of objects of similar magnitudes as our targets. We cor- 
rect for the intrinsic dispersion based on the Besancon 
model simulation, although we estimate that measure- 
ment errors are the main contributor to the total dis- 
persion. These two arguments suggest that our error 
estimates are correct within a few mas/yr. 

In that case, the proper motion of CAHA Tau 1 is 
incompatible with a giant star proper motion at the l.Scr 
level; CAHA Tau 2 at the 1.9ct level. We note that the fit 
of CAHA Tau 2 proper motion along the right ascension 
is not good, and a value larger than Taurus members' 
proper motions cannot be confidently excluded. 

3.4.2. CAHA Tau 3, 4 and 5 

Unfortunately the situation for CAHA Tau 3, CAHA 
Tau 4 and CAHA Tau 5 is different. The data set for 
these targets is more limited than for CAHA Tau 1 and 
CAHA Tau 2, in particular no UKIDSS data are avail- 
able. Also the positional dispersions of reference stars is 
larger, pointing to a poorer registration. Based on the 
current data we find that the proper motions we mea- 
sure for CAHA Tau 3 do not correspond to the typical 
proper motions of known Taurus members (see Table A6 
and Figure AlO). However, the current error bars are 
large and our measurements are not stable against the 
registering parameters' fine tuning, so that further ob- 
servations are required to ensure a robust measurement. 
The proper motions of CAHA Tau 4 and CAHA Tau 5 
are in agreement with a Taurus membership but have 
little contamination rejection power. 

4. DISCUSSION 
4.1. CAHA Tau 1 - the first L2 dwarf in Taurus 

From the results shown above we believe that CAHA 
Tau 1 is a young, low-mass object in Taurus with spec- 
tral type L2±0.5 and thus so far the object with the 
latest spectral type in Taurus (Luhman et al. 2009a) . To 
estimate its effective temperature we used the relation 
between spectral type and Tg// for low-mass field ob- 
jects found by Golimowski et al. (2004). Based on the 
same reference we also derived the bolometric correction 
in the K-h&nd. All values are summarized in Table A7. 
The errors are based on the uncertainties in the spectral 
type and the intrinsic errors of the conversions found by 
Golimowski et al. (2004). 

In Figure A7 we plot an HR-diagram for CAHA Tau 1 
and known low-mass Taurus objects with spectral types 
>M9.0. The luminosity of CAHA Tau 1 is derived from 
the bolometric correction BCk (Table A7) and assuming 
a distance to Taurus of 140 pc and Av='2.Z. The refer- 
ence objects were collected from Luhman et al. (2009b, a) 
and Guieu et al. (2006). Overplotted in the diagram are 
isochrones (1, 5, 10, 50, and 100 Myr) and mass tracks 
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(5, 10, 15, 30, and 50 Mjup) from Chabricr ct al. (2000). 
Three points arc noteworthy: (1) For most objects the 
derived properties are in good agreement with young (< 
10 Myr), low mass (< 30 Mj„p) objects. (2) CAHA 
Tau 2 (number 1 in the plot) falls nicely between the 
1 and 5 Myr isochrones and, depending on its actual 
age, its observed properties translate into a most likely 
mass between 5 and 15 Mj„p based on theoretical mod- 
els (Table A7). It is thus the coolest Taurus member 
known today and only the second object in the region 
with a mass most likely below the Deuterium burning 
limit making it a free floating planetary mass object. 
(3) The other object that was recently identified as a 
Taurus member of planetary mass with only 4-7 Mjup 
(Luhman et al. 2009a, 2MASS 0419-1-2712, number 2a 
in the plot) lies below the 100 Myr isochrone and close 
to the 50 M,7up mass track. This behavior was already 
discussed in the discovery paper where the authors con- 
cluded that the old isochronical age for the object is likely 
a reflection of errors in the adopted temperature scale 
and evolutionary tracks. Similarly, the coolest member 
of the Chameleon I SF region (spectral LO) appears also 
too old for its luminosity (Luhman et al. 2008) which is 
comparable to that of 2MASS 0419+2712. Disregarding 
the derived temperature and assuming an age of 1-3 Myr 
the derived luminosity of 2MASS 0419-1-2712 yields then 
a mass of 4-7 M,7„p as mentioned above. While there are 
certainly notable uncertainties in the temperature scales 
and evolutionary models, it is interesting to note that 
all other late type Taurus members shown in Figure A7 
(inch CAHA Tau 1) populate as expected the young 
areas of the diagram. Thus, it is worth investigating 
possible astrophysical reasons why 2MASS 0419-1-2712, 
if it is indeed a Taurus member, should appear under- 
luminous. As a flrst test, we use the same temperature 
scale and bolometric correction that we used for CAHA 
Tau 1 . This leads indeed to a significant shift in the HR- 
diagram (number 2b): Instead the of original values of 
Te//=2200 K and log L/Lq w3.52, we find re//=2300 
K and log L/L© «3.26, demonstrating the uncertain- 
ties across the diff'erent temperature scales. However, 
2MASS 04194-2712 appears stiU too low in the diagram 
compared to the other sources. As a second source of 
uncertainty, we speculate that part of this might result 
from underestimating the extinction towards the source. 
Even if no disk, and in particular no edge-on disk creat- 
ing significant extinction, appears to be present (Luhman 
et al. 2009a), the assumed value of Aj^O mag, hence no 
extinction, for 2MASS 0419+2712 might be too low. We 
note that also in the case of the LO dwarf in Chameleon I 
(Luhman et al. 2008) no extinction was assumed. Finally, 
we mention that it could still be possible that 2MASS 
0419+2712 is indeed older than typical Taurus members 
or it is not associated with the main Taurus population 
and further away than the assumed 140 pc. In all these 
case the mass of the object would be higher than the 
current estimates of a few Jupiter masses. More spec- 
troscopic data and a proper motion analysis might yield 
additional insight into this puzzle. 

As CAHA Tau 1 fits nicely on the HR diagram we 
compare its broad band SED with simulations. In Fig- 
ure A8 we plot the dereddened photometry together with 
a theoretical model SED for a 3 Myr object of 10 Jupiter 
masses with Te//=2040 K (Burrows et al. 2003, 2006). 



The fiuxes were normalized to unity in the _ff-band. The 
observed photometry was corrected for extinction using 
the value for Ay given in Table A5 and the extinction 
curve from Mathis (1990) with Rv=3.1. The agreement 
between the model and the observations is very good and 
only the observed K-band and 4.5 fim fluxes appear a lit- 
tle too high compared to the simulations. At least for the 
flux at 4.5 ^m it might well be that circum(sub-)stellar 
material creates excess emission above the photosphere. 
Zapatero Osorio et al. (2007) found also excess emission 
for their young, very low-mass objects but mostly long- 
ward of 5 /im. It would be nice to check whether the 
apparent emission excess we see in CAHA Tau 1 extends 
also to longer wavelengths. 

4.2. CAHA Tau 2 - 5 - Embedded carbon stars or 

(super-) giants? 

While for CAHA Tau 1 there are sufficient indica- 
tions that it is a young Taurus member, it has not yet 

been convincingly shown that our candidates 2-5 are 
not background carbon stars or late-type (super-)giants. 
And eventually only additional spectra will help to settle 
this open question. However, to elaborate a little further 
on this problem, we plot in Figure A9 a color-magnitude 
diagram with our candidates and the two groups of car- 
bon stars we already mentioned above (Tanaka et al. 
2007; Liebcrt et al. 2000). For completeness also CAHA 
Tau 1 and CAHA Tau 6 (alias FU Tau B) are shown 
in the plot. One finds that the very distant and deeply 
embedded objects from Liebcrt et al. (2000) tend to be 
much redder than our candidates while the nearby car- 
bon stars from Tanaka ct al. (2007) have a very similar 
J — Kg color like the bulk of our objects. Even more 
important, however, is the fact that all of our candi- 
dates appear fainter than even the most distant objects 
from Liebcrt et al. (2000) which are already assumed 
to lie in or even beyond the Galactic halo. Hence, if 
our candidates were distant objects then (1) they should 
appear redder than they are as the Taurus molecular 
cloud should cause significant additional extinction and 
(2) they should be extremely far away. In Table A8 we 
list the typical values for absolute K-band magnitudes of 
carbon stars and M-type giants. Based on these figures 
and the apparent magnitudes of our remaining 2 least 
luminous candidates (CAHA Tau 2 and 3; Table A5) we 
derive a range for the objects' potential distance moduli 
{m,K — Mk)- We focus again on these objects as they 
are intrinsically more interesting due to their faintness. 
Using 

d = lQ0.2-{mK-MK) + l (^i^ 

we can estimate the physical distances d [in pc] for our 
candidates if they were indeed carbon stars or M-giants. 
It shows that the minimum distances for these objects 
would be around -^155 kpc for carbon stars and ~135 
kpc for late type M giants (Table A8). If the assumed ex- 
tinction towards the objects (Table A5) is not underesti- 
mated the derived distances lie way beyond the currently 
known extent of our Galaxy, and to our knowledge car- 
bon stars have only been found at distances up to ~130 
kpc (Mauron et al. 2007, and references therein). Only 
the minimum distance value for early M-typc giants ^60 
kpc appears still reasonable, but in Figure A5 we showed 
that at least CAHA Tau 2 is not an early M-type object. 
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The density law of halo (giant) stars at those large 
Galacto-centrie distances is largely unconstrained, so 
that it is difficult to estimate even roughly the number 
of such stars in our total data sample. Extrapolating 
the halo giant stellar count of Majewski et al. (2003) as- 
suming a smooth halo and a surface density law, 
we expect much less than one object in our survey area. 
However, it is possible that the halo be highly struc- 
tured at these distances, so that the counts would be 
dominated by (unknown) substructures. In particular, 
we note that Taurus lies on the Sagittarius stream or- 
bital plane, on its Southern arc. However, no Sagittarius 
tails are known past « 40kpc, or predicted past 100 kpc 
(Law et al. 2005). 

We think that Figure A9 and Table A8 provide ad- 
ditional arguments that most of our objects can be re- 
garded as potential low-mass candidates. Especially the 
least luminous objects CAHA Tau 2 and 3 would be ex- 
tremely distant if they were background objects. The 
only exception appears to be CAHA Tau 4. Already its 
photometric properties discussed in section 3.2 showed 
that this object might be of different nature than the 
other candidates. Figure A9 further strengthens this hy- 
pothesis as CAHA Tau 4 could potentially be a (very) 
distant carbon star. 

4.3. Comparison to theoretical models 

Even though only spectra can eventually prove the 
young and low-mass nature of CAHA Tau 2 - 5 we pro- 
vide in the following a first c;omparison of their photo- 
metric properties to theoretical models. The broadband 
SED of CAHA Tau 1 was already discussed above (Fig- 
ure AS) and that of CAHA Tau 6 (alias FU Tau B) was 
analyzed in Luhman et al. (2009b) . For the other objects 
we derived estimates for effective temperatures, potential 
masses and surface gravities from the models of Burrows 
et al. (2003, 2006) (for an assumed age of 3 Myr) and 
Chabrier et al. (2000) (for an assumed age of 1 Myr). 
To derive these parameters we fitted the extinction cor- 
rected photometric values for the /, R and J band (Ta- 
ble A5) to model grids covering a mass regime between 
5 and 70 Mj„p. The I, R and J band were chosen for 
the following reasons: (1) for these bands suitable model 
grids covering the mass regime were accessible for both 
ages; (2) focusing on the shorter wavelength regime of 
our observations reduces the potential contamination by 
circumstellar material altering the NIR fluxes; (3) the J- 
band photometry has a smaller intrinsic error than the 
optical wavelength bands and add an additional color to 
our fits^^. For the extinction correction we used the in- 
dividual values for Ay given in Table A5 and applied the 
extinction law from Mathis (1990). We furthermore as- 
sumed a distance of 140 pc (corresponding to a distance 
modulus of 5.73 mag) and adapted the photometric val- 
ues of the theoretical models (originally given in absolute 
magnitudes) accordingly. 

Table A9 summarizes the fit results assuming that our 
objects are indeed young, low-mass objects. Changing 
the distance to 160 pc (performed as an additional test 
case) did only alter the results for one object, CAHA Tau 

The derived best fit masses did not change for most objects 
when we used only the / and R band. Only for CAHA Tau 3 we 
found slightly lower masses in this case. 



2, in the 1 Myr case. Here, a slightly higher mass was 
found compared to the value given in Table A9. Over- 
all, this fitting exercise shows that, depending on the 
age, the first three objects would be among the least 
massive and coolest Taurus objects known today in case 
they were Taurus objects. In particular, CAHA Tau 2 
cotild lie close to the deuterium burning limit and thus 
on the edge to planetary mass objects. However, as al- 
ready mentioned in section 3.3, only spectra will yield 
the ultimate answer to spectral type and thus mass of 
our objects. 

In Figure A12 we compare all observed, extinction cor- 
rected magnitudes to the theoretical broadband SEDs of 
the best fit models given in Table A9. It shows that the 
magnitudes for CAHA Tau 2 and 4 at the shorter wave- 
length end are already in general agreement with the 
models while the restilts for CAHA 3 and 5 show notable 
deviations. Even if we imply that the objects are young, 
low-mass objects, as done in this Figure, then there are 
still two major sources of uncertainty: (1) the individ- 
ual extinction can differ from the value derived from the 
extinction maps, and (2) the distance, and thus the re- 
quired distance modulus, of the objects may vary. The 
plot for CAHA Tau 5, and partly also that for CAHA Tau 

3, yield hints that the extinction might have been over- 
estimated as the slope between the three optical SDSS 
bands is not as steep as required by any of the mod- 
els. However, we emphasize again that only with the 
help of additional spectra the currently existing degen- 
eracy between spectral type (and thus mass), Ay and 
distance can be resolved and we refrain from further re- 
finements of the exercise presented here. Also, the ques- 
tion of NIR/MIR excess emission due to a disk can only 
be answered once the spectral type and Ay is confirmed 
for each object. For instance, while CAHA Tau 4 shows 
already with the current assumptions NIR excess emis- 
sion, CAHA Tau 2 shows only excess emission compared 
to the 3 Myr models. 

5. CONCLUSIONS 

Based on a deep NIR survey in Taurus we identified 6 

new candidates for low-mass objects in this young star- 
forming region. One object, CAHA Tau 6 (now FU Tau 
B), was independently confirmed as young, low-mass BD 
of spectral type M9.25 by Luhman et al. (2009b). For a 
second object CAHA Tau 1, we have spectroscopic ev- 
idence that it is indeed a young, very low-mass objects 
and we assigned a spectral of L2±0.5, making it the ob- 
ject with the latest spectral type in Taurus. For the 
remaining 4 objects we thus far lack clear spectroscopic 
confirmation for their young, low-mass nature and we 
can only consider them as candidates at this point in 
time. However, the identification of CAHA Tau 1 and 
CAHA Tau 6 with spectral types L2 and M9.25, respec- 
tively, as well as the re-identification of 10 additional 
known low-mass Taurus members during the selection 
process, gives us confidence for the low-mass nature of 
the other objects. If they turn out to be indeed young 
Taurus sources then a comparison to theoretical models 
suggests that at least a second object (CAHA Tau 2) 
could be among the least massive known objects in this 
young star-forming region known today. In this case, its 
mass would be close to or below the deuterium burning 
limit of ~13 Mjup and its effective temperature would 



Very low-mass brown dwarfs and planetary-mass objects in Taurus 



9 



lie below ~2300 K making it also an L-type object. The 
findings in detail: 

• CAHA Tau 1: The photometric properties, its 
proper motion and also the shape and spectro- 
scopic features in the J-band spectrum make it 
very likely that CAHA Tau 1 is the first Taurus 
member with a spectral type as late as L2±0.5 and 
an effective temperature of ^^2100 K. For an age 
between 1 and 5 Myr, as derived from its position 
in a HR-diagram, this temperature corresponds to 
5-15 Mjup and it might well be that CAHA Tail f 
is a free floating planetary mass object. Its broad 
band SED is consistent with model predictions for 
a 3 Myr-old object of 10 Jupiter masses with slight 
excess emission at 4.5 /xm possibly arising from a 
circum(sub-)stellar disk. 

• CAHA Tau 2: The optical and NIR photometric 
properties support the young, low-mass hypothe- 
sis for this object. Further evidence comes from 
its proper motion. A first J-band spectrum did 
not unambiguously prove the young, low-mass na- 
ture as the long wavelength end resembles that of 
evolved stars or carbon stars. However, a compari- 
son of the photometric properties makes this rather 
unlikely. Assuming that CAHA Tau 2 is young, 
low-mass objects then theoretical models suggest 
that, depending on the age, it lies close to or be- 
low the deuterium burning limit making it another 
candidate for being among the least massive and 
coolest known Taurus members today. 

• CAHA Tau 3: For CAHA Tau 3 we lack spectro- 
scopic information but its photometric properties 

are in good agreement with those of yoiuig, very 
low-mass objects. However, currently the proper 
motion analysis does not support a Taurus mem- 
bership. Additional data is needed to reduce the 
large error bars and get a more robust assessment 
of the proper motion. Assuming a young, low-mass 
nature then theoretical models suggest that CAHA 
Tau 3 lies somewhere at the border between M- and 
L-type objects with ~15 M_Jup. In this case, the 
currently applied value for the extinction is prob- 
ably overestimated and needs to be revised which 
would again reduce the estimated mass. 

• CAHA Tau 4 and 5: The optical and NIR data for 
CAHA Tau 5 are in good agreement with young 
BD while CAHA Tau 4, the brightest object in our 
sample, lies a little bit off in the presented color- 
color and color-magnitude plots. Based on its pho- 
tometric properties CAHA Tau 4 could also be a 



(very) distant, evolved object. A fit to theoretical 
models for young low-mass objects shows that if 
CAHA Tau 4 was a young BD then it would prob- 
ably have significant NIR excess emission. Fitting 
models to the photometric values of CAHA Tau 
5 reveals that the currently applied value for the 
extinction (and/or possibly also the applied dis- 
tance modulus) need to be refined once spectro- 
scopic data is available. Our PM analysis does not 
put conclusive constraints on the Taurus member- 
ship for both sources. 
• CAHA Tau 6: This object was identified as FU Tau 
B, a young, M9.25 object in a BD binary system 
by Luhman et al. (2009b) during the referee process 
of this paper. We publish here the first J, H and 
Kg band data for this source complementing the 
broadband SED from Luhman et al. (2009b). 

So far we have not yet made an attempt to estimate the 
effect of our findings on the Taurus IMF. Once we know 
better which of the above mentioned candidates are in- 
deed young, low-mass sources and how many additional 
candidates we find in the remaining and not yet here in- 
cluded fields of our deep NIR survey, it is certainly worth 
analyzing whether our results significantly influence the 
low-mass end of the IMF and help to fill-up the apparent 
deficiency of very low-mass objects compared to late K 
and early M type Taurus members. 
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APPENDIX 
TRANSFORMATION EQUATIONS 
2MASS - GIT 

(see, http: / /www. ipac.caltech.edu/2mass/relcascs/allsky/doc /sec6_4b.html) 

(i^«)2MASS = KciT - (0.019 ± 0.004) + (0.001 ± 0.005)( J - K)cn 
(J-H) 

2MASS = (1.087 ± 0.013)( J - iJ)ciT - (0.047 ± 0.007) 
(J - ^a)2MASS = (1-068 ± 0.009)(J - K)ciT - (0.020 ± 0.007) 
{H - /<:«)2MASS = (1.000 ± 0.023) (if - K)ciT + (0.034 ± 0.006) 
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2MASS - BESSEL & BRETT 
(see, http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4b.html) 

(^«)2MASS = ^BB - (0.039 ± 0.007) + (0.001 ± 0.005)(J - K)^^ 
{J - ff)2MAss = (0.990 ± 0.012)(J - ff)BB - (0.049 ± 0.007) 
{J - Ks) 2MASS = (0.983 ± 0.008) (J - i^)BB - (0.018 ± 0.007) 
{H - X.)2MAss = (0.971 ± 0.022)(i/ - K)bb + (0.034 ± 0.006) 

SDSS (UGRIZ; - JOHNSON-COUSINS (UBVRi; 

(see, Jordi et al. 2006) 

r-z = (1.584 ± 0.008) (i? -I)- (0.386 ± 0.005) 
i-I= (0.247 ± 0.003) (i? - /) + (0.329 ± 0.002) 
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Quanz et al. 



TABLE Al 

Estimated masses of objects with an effective 
temperature of ~2500 k for different ages based on 
theoretical models. masses are given in jupiter 

MASSES. 



Model 


1 Myr 


2 Myr 


3 Myr 


5 Myr 


BarafFe et al. (2003) 
Chabrier et al. (2000) 
Burrows et al. (2003, 2006) 


~20 Mjup 
~17 Mj„p 


~16 Mj^p 


~17 Mj^p 


~23 Mj„p 
~20 Mjup 
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TABLE A2 

Masses and effective temperatures for young, low-mass objects of different ages 

BASED ON THEORETICAL ^[OnELS. MASSES ARE On'EK IX JrriTER MASSES. 



Model 1 Myr 2 Myr 3 Myr 5 Myr 



Baraffe et al. (2003) 5 Mj„p / 1900 K - - 5 Mj„p / 1410 K 

10 Mj„p / 2250 K - - 10 Mj„p / 1970 K 

15 Mjup I 2400 K - - 15 Mj^p / 2260 K 

Chabrier et al. (2000) 5 Mj^p / 1900 K - - 5 Mj„p / 1410 K 

10 Mj„p / 2240 K - - 10 M^juv I 1870 K 

15 Mj„p / 2450 K - - 15 M.j^^ / 2250 K 

Burrows et al. (2003, 2006) - 5 Mj„p / 1620 K 5 M / 1500 K 

10 Mj„p / 2130 K 10 Mj„p / 2040 K 

15 Mjup I 2450 K 15 Mj„p / 2390 K 
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TABLE A3 

Journal of the Omega2000 observations in Taurus. Column four shows the number of 1-minute exposures per field and filter 

that were finally co-added to create the final images. the seeing (column 5) was estimated during the data reduction in 
the respective nir filters. the sixth column shows the mean transformation error for each field and filter resulting from 

the photometric calibration with 2mass data. 



Field name 


RA (J2000) 


DEC (J2000) # of exposures 


Seeing ["] 


Calibr. err. [mag] 


Obs. date 


SDSS 






J/H/Ks 


J/H/Ks 


J/H/Ks 




coverage 



tc2 


04hl0m53, 


.64s 


+25°09'28.3" 


28 


/ 


28 


/ 


30 


1 


.3 


/ 


1. 


.0 


/ 


1, 


.0 


0, 


.04 


/ 


0, 


.05 


/ 


0, 


.05 


Jan. 


05 


50% 


tc5.6.7.8 


04hl8m00, 


.87s 


+28°12'28.9" 


30 


/ 


28 


/ 


30 


1 


.0 


/ 


1, 


.0 


/ 


1, 


.5 


0, 


.07 


/ 


0, 


.05 


/ 


0, 


.05 


Jan. 


05 


100% 


tc7.11 


04hl8ml9, 


.56s 


+28°23'59.2" 


30 


/ 


30 


/ 


30 


1 


A 


/ 


1. 


.0 


/ 


1, 


.3 


0, 


.09 


/ 


0, 


.07 


/ 


0, 


.07 


Oct. 


04 


100% 


tcl2a_12b 


04hl9m03, 


.83s 


+27°19'27.4" 


30 


/ 


30 


/ 


21 


1 


.6 


/ 


1. 


.7 


/ 


1, 


.7 


0, 


.05 


/ 


0, 


.03 


/ 


0, 


.06 


Oct. 


04 


92% 


tcl3a_13b_14 


04hl9m44, 


.22s 


+27°09'52.2" 


30 


/ 


27 


/ 


29 


1 


.2 


/ 


1, 


.0 


/ 


1, 


.2 


0, 


.06 


/ 


0, 


.07 


/ 


0, 


.05 


Jan. 


05 


100% 


tcl5_16a 


04h20m35, 


.70s 


+27°05'37.2" 


30 


/ 


30 


/ 


30 


1 


.4 


/ 


1, 


.3 


/ 


1 


.3 


0, 


.08 


/ 


0, 


.08 


/ 


0, 


.08 


Dec. 


05 


100% 


tcl6a.l6b.l7 


04h21m22, 


.49s 


+26°59'08.0" 


30 


/ 


30 


/ 


30 


1 


.6 


/ 


1, 


.4 


/ 


1 


.3 


0, 


.09 


/ 


0, 


.05 


/ 


0, 


.07 


Dec. 


05 


100% 


tcl9 


04h23m36, 


.53s 


+25°05'55.3" 


30 


/ 


30 


/ 


30 


1 


.3 


/ 


1, 


.3 


/ 


1, 


.2 


0, 


.06 


/ 


0, 


.04 


/ 


0, 


.05 


Dec. 


05 


100% 


tc20.21_22 


04h23m55, 


.75s 


+26°36'38.4" 


30 


/ 


30 


/ 


30 


1 


.8 


/ 


1, 


.7 


/ 


1, 


.2 


0, 


.05 


/ 


0, 


.05 


/ 


0, 


.06 


Oct. 


04 


100% 


tc23.24 


04h26m33, 


.31s 


+24°38'47.4" 


30 


/ 


30 


/ 


29 


1 


.4 


/ 


1, 


.5 


/ 


1, 


.3 


0, 


.06 


/ 


0, 


.04 


/ 


0, 


.05 


Oct. 


04 


100% 


tc25.26a_26b 


04h27m53, 


.53s 


+26°19'01.0" 


30 


/ 


30 


/ 


30 


1 


.2 


/ 


1, 


.6 


/ 


1 


.2 


0, 


.04 


/ 


0, 


.06 


/ 


0, 


.05 


Dec. 


05 


100% 


tc28.29 


04h29m40, 


.00s 


+24°29'34.1" 


30 


/ 


30 


/ 


28 


1 


.9 


/ 


2, 


.4 


/ 


2, 


.4 


0, 


.06 


/ 


0, 


.04 


/ 


0, 


.06 


Oct. 


04 


100% 


tc35_36a_36b 


04h35m44 


.37s 


+24°09'11.8" 


29 


/ 


30 


/ 


30 


1 


.2 


/ 


1, 


.1 


/ 


1 


.0 


0, 


.06 


/ 


0, 


.06 


/ 


0, 


.06 


Jan. 


05 


77% 


tc37_38.39.41 


04h39m28, 


.62s 


+25°47'29.6" 


29 


/ 


29 


/ 


30 


1 


.0 


/ 


1, 


.0 


/ 


1, 


.1 


0, 


.09 


/ 


0, 


.07 


/ 


0, 


.05 


Jan. 


05 


100% 


tc42a-42b.42c 


04h40m36, 


.19s 


-|-25°29'53.4" 


29 


/ 


30 


/ 


30 





.9 


/ 


1, 


.0 


/ 


0, 


.9 


0, 


.07 


/ 


0, 


.05 


/ 


0, 


.05 


Jan. 


05 


100% 


tc43.44 


04h41m28, 


.61s 


+25° 54' 15.5" 


30 


/ 


29 


/ 


30 


1 


.0 


/ 


1, 


.1 


/ 


1, 


.0 


0, 


.06 


/ 


0, 


.05 


/ 


0, 


.10 


Jan. 


05 


100% 


tc50 


04h41m30 


.00s 


+25°42'30.0" 


30 


/ 


30 


/ 


30 


1 


.2 


/ 


1. 


.2 


/ 


1, 


.1 


0. 


.07 


/ 


0. 


.06 


/ 


0. 


.06 


Dec. 05/ Jan. 06 


100% 



tcla_lb 
tc3.4 

tc9.10a.10b 
tcl8a.l8b 
tc27 
tc30 

tc31.32.33a.33b 
tc34 



04h04m46.45s 
04hl3m52.10s 
04hl8m25.60s 
04h21m58.98s 
04h28m39.26s 
04h31m41.90s 
04h32ml9.18s 
04h33m24.28s 



+26° 19' 27.7" 
+28° 15' 45.8" 
+27°30'01.3" 
+15°30'02.4" 
+26°51'39.2" 
+18°09'50.7" 
+24°28'00.6" 
+22°42'15.9" 



30 / 43 / 30 
30 / 30 / 30 
30 / 30 / - 
30 / 29 / 28 
30 / 30 / 30 
30 / 30 / 30 
30 / 30 / - 
28 / 29 / 27 



1.7 / 1.4 / 1.3 
1.0 / 1.1 / 1.0 
2.3 / 1.8 / - 

1.1 / 1.2 / 1.3 

1.5 / 1.3 / 1.3 

2.6 / 1.9 / 2.4 
1.2 / 1.3 / - 

1.2 / 1.2 / 1.2 



0.11 / 0.05 / 0.06 
0.05 / 0.04 / 0.03 
0.08 / 0.04 / - 
0.04 / 0.04 / 0.06 
0.07 / 0.07 / 0.06 
0.07 / 0.06 / 0.05 

0.06 / 0.06 / - 
0.06 / 0.06 / 0.07 



Dec. 05 
Jan. 05 
Dec. 05 
Jan. /Feb. 05 
Dec. 05 
Oct. 04 
Dec. 05 
Jan. /Feb. 05 
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TABLE A4 

Properties of selected but already known low-mass objects in Taurus. For each object magnitudes from our 
observations as well as from the 2mass catalogue (second line of magnitudes) are given for comparison. the 
errors in our nir magnitudes reflect uncertainties in the fitting of the individual psfs as well as systematic 

ERRORS RESULTING FROM THE PHOTOMETRIC CALIBRATION WITH THE 2MASS REFERENCE SOURCES. REFERENCES: (a) LuHMAN 

ET AL. (2006a), (b) Luhman (2006), (c) GuiEU ET AL. (2006), (d) White & Basri (2003), (e) Luhman et al. [1998), (f) 

MOHANTY ET AL. (2005), (c) HARRIS ET AL. (1988), (h) ItOH ET AL. (1999), (l) ItOH ET AL. (2002), (j) BrICENO ET AL. 

(2002), (k) 2MASS Point Source Catalogue (Cutri et al. 2003), (l) Luhman et al. (2009b) 



Object 



RA (J2000) DEC (J2000) 



H 



Ks 



Spec. Type 



Ref. 



V410 X-ray 3 04hl8m08.0s -|-28°26'03.9" 

IRAS 04154-1-2823 04hl8m32.1s -|-28°31'15.5" 

2MASS J04242090-I-263051 04h24m20.9s -f26°30'51.0" 

2MASS J04263055+2443558 04h26m30.6s -f24°43'55.8" 

CFHT-BD-Tau 20 04h29m59.5s +24°33'07.7" 

ITG33 04h41m08.3s -f25°56'07.3" 

CFHT-BD-Tau 19 04h21m08.0s +27°02'20.1" 

2MASS J04215450-f2652315 04h21m54.5s -t-26°52'32.0" 

FU Tau A 04h23m35.4s -f25°03'02.5" 

KPNO-Tau 4 04h27m28.0s -|-26° 12' 04.7" 



11.60±0.10 
11.54±0.02 
14.38±0.10 
15.19±0.05 
13.40±0.06 
13.49±0.02 
14.74±0.07 
14.67±0.03 
11.77±0.07 
11.68±0.02 
14.01±0.07 
13.74it0.03 
12.00±0.10 
13.86±0.02 
15.53±0.10 
15.54±0.04 
10.77±0.07 
10.78±0.03 
15.06±0.05 
15.00±0.04 



10.91±0.09 
10.82±0.02 
12.27±0.10 
12.37±0.02 
12.82±0.06 
12.81±0.02 
13.95±0.05 
13.94±0.04 
10.89±0.05 
10.54±0.03 
12.51±0.08 
12.15±0.02 
10.33±0.09 
12.06±0.03 
14.46±0.06 
14.50±0.04 
9.98±0.06 
9.95±0.03 
13.97±0.07 
14.03±0.04 



10.52±0.09 
10.45±0.02 
9.59±0.08 
10.27±0.02 
12.48±0.07 
12.43±0.02 
13.40±0.06 
13.40±0.03 
10.14±0.07 
9.81±0.02 
11.62±0.13 
11.09±0.02 
9.19±0.12 
10.54±0.02 
13.81±0.09 
13.90±0.04 
9.45±0.07 
9.32±0.02 
13.29±0.06 
13.28±0.03 



M6-M6.5 

M6.5 
M8.75 

M5.5 
K7-M3 
M5.25 

M8.5 
M7.25 

M9.5 



(d),(e),(f) 

(k) 

(g) 

(k) 

(a) 

(k) 

(b) 

(k) 

(c) 

(k) 
{h),(i) 

(k) 

(c) 

(k) 

(b) 

(k) 
{g)/{l) 

(k) 

(j) 

(k) 
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TABLE A5 

Coordinates and photometry of the newly identified candidates for young (very) low-mass objects in 
Taurus as well as rough estimates for their optical extinction. Optical magnitudes are given in the SDSS 
photometric system (r, i, z) and johnson-cousin system {r,i), the nir magnitudes (j, h, kg) are in the 

2MASS SYSTEM, AND THE LONGEST WAVELENGTHS GIVEN HERE ARE SPITZER/IRAC MAGNITUDES. THE ERRORS IN THE R 
AND / BAND ARE DERIVED FROM THE TRANSFORMATION EQUATION (SEE APPENDIX) AND THE INDIVIDUAL ERRORS IN THE 
SDSS MAGNITUDES. THE NIR ERRORS ACCOUNT FOR TRANSFORMATION ERRORS BETWEEN THE OmEGA2000 FILTER 
SYSTEM AND THE 2MASS PHOTOMETRIC SYSTEM AS WELL AS FOR INDIVIDUAL FITTING ERRORS DURING THE PSF 
PHOTOMETRY. THE ERRORS IN THE IRAC BANDS REFLECT THE UNCERTAINTIES IN THE APERTURE PHOTOMETRY. 



Object 



RA (,12000) DEC (J2000) 



R 



CAHA Tau 1 
CAllA Tau 2 
CAHA Tau 3 
CAHA Tau 4 
CAHA Tau 5 



04h36m09.0s 
04h39m47.2s 
04h40m39.3s 
04h40m57.4s 
04h41m31.6s 



+24°08'21.2" 
+25°53'32.2" 
+25°23'02.0" 
+25°50'07.9" 
-|-25°49'32.2" 



23.42±0.41 
22.65±0.22 
23.60±0.48 
23.41±0.41 
20.17±0.03 



21.09±0.08 
20.17±0.04 
21.47±0.11 
20.74±0.07 
18.68±0.01 



19.34±0.07 
18.66±0.04 
19.73±0.12 
18.81±0.05 
17.47±0.02 



22.89±0.38 
21.92±0.22 
23.16±0.57 
22.78±0.41 
19.82±0.05 



20.00±0.35 
19.16±0.17 
20.47±0.47 
19.63±0.31 
17.87±0.04 



CAHA Tau 6* 04h23m35.8s +25°02'59.4" 22.25±0.14 20.19±0.05 17.80±0.05 22.16±0.22 19.11±0.20 



H 



3.6 jXTa 



4.5 nm. 



5.8 ixm. 



CAHA Tau 1 
CAHA Tau 2 
CAHA Tau 3 
CAHA Tau 4 
CAHA Tau 5 
CAHA Tau 6" 



16.92±0.07 
16.33±0.12 
17.46±0.10 
15.75±0.07 
14.90±0.14 
15.10±0.08 



15.77±0.09 
15.20±0.10 
16.13±0.08 
14.12±0.08 
14.39±0.08 
14.12±0.05 



15.11±0.09 
14.55±0.08 
15.12±0.10 
13.04±0.17 
13.90±0.08 
13.43±0.07 



14.62±0.08 
14.21±0.07'^ 
14.84±0.08 
12.90±0.06 
13.65±0.06 
12.54±0.1= 



14.43±0.09 
14.14±0.08'^ 
14.83d=0.11 
12.82±0.06 
13.53±0.07 
11.90±0.1= 



12.73±0.16 
13.46±0.26 
11.46±0.1'= 



~ 2.3 

~ 0.9 
~ 4.1 
~ 4.7 
~ 4.0 
' 4.9/< 1.0== 



^Identified as FU Tau B by Luhman et al. (2009b) during the revision process of this paper. 

''Object lies in so-called "pull-down" column from a bright star in the IRAC images. IRAC fluxes axe probably under- 
estimated. 

"^Figures taken from Luhman et al. (2009b). 

■^Values are estimates based on the Taurus extinction map by Dobashi et al. (2005) which is created from DSSl optical 
data with a resolution of ~ 6'. 
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TABLE A6 

Proper motions and time baseline of our candidates. 



IMa [mas/yr] /x^ [mas/yr] time baseline [d] 



CAHA Tau 1 +12 ± 20 -26 ± 15 1373 

CAHA Tau 2 +10 ± 15 -28 ± 10 1566 

CAHA Tau 3 -42 ±40 +7 ± 20 1355 

CAHA Tau 4 +29 ±20 -17 ± 30 2270 

CAHA Tau 5 -8 ± 16 -10 ±16 1514 
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TABLE A7 

Properties of CAHA Tau 1 based on its J-band spectrum. 



Spectral type Tj.// [K]'^ BCk"" log L/L©'' Age [Myr] Mass [Mj„p]<= 



L2±0.5 2080±140 3.3±0.2 3.1±0.1 1/3/5 7+1 / WZl / 12; 



''Based on the relations between, SpT, T^// and bolometric correction in the K band 
(BCx) found by Golimowski et al. (2004). 

''Derived using: L = 3.02 ■ 10^* X 10-° '* *^f'o( [w] (Lang 1980). 

"^Most likely mass range for an age of 1, 3 and 5 Myr based on models of Chabrier 
et al. (2000) and Burrows et al. (2003, 2006). 
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TABLE A8 

Typical absolute K-band magnitudes of carbon stars and M-type 
giants and corresponding minimum and maximum distance moduli and 

distances if our REMAINING TWO LEAST LUMINOUS CANDIDATES (CAHA TAU 2 
AND 3) WERE SUCH OBJECTS. CAHA TAU 2 SERVES AS BASIS FOR THE MINIMUM 

ESTIMATE {rriK = 14.45), CAHA Tau 3 for the maximum (niK = 14.67)''. 
References: (a) Liebert et al. (2000); Claussen et al. (1987), (b) 
SIMBAD Astronomical Database (2MASS catalogue, Hipparcos 
catalogue): HD 196610, HD 18191, HD 207076, HD 108849; (c) SIMBAD 
Astronomical Database (2MASS catalogue, Hipparcos catalogue): HD 
213893, HD 204724, HD 120052, HD 219734 



Object type 


Mk [mag] 
min / max 


Reference 


txik-Mk [mag]" 
min / max 


Distance [kpc] 
min / max 


Carbon stars 
M6 HI - M8 HI 
MO HI - M2.5 HI 


-6.5 / -8.5 
-6.2 / -7.4 
-4.2 / -6.0 


(a) 
(b) 
(c) 


20.95 / 23.17 
20.65 / 22.07 
18.65 / 20.67 


155 / 431 
135 / 259 
54 / 136 



''The apparent X-band magnitudes of CAHA Tau 2 and 3 were corrected for 
extinction using the optical extinction values given in Table A5 and assuming 
Ak/Av s» 0.11 (Mathis 1990). 
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TABLE A9 

Best fit results for mass, effective temperature and surface gravity of 

OUR candidates for 1 AND 3 MyR. THE FIGURES FOR THE EXTINCTION ARE THE 

SAME AS IN Table A5. The left figure in each column reflects the best 

FIT TO the 1 MyR models BY CHABRIER ET AL. (2000), THE RIGHT FIGURE IN 
each COLUMN REFLECTS THE BEST FIT TO THE 3 MYR MODELS BY BURROWS 

ET AL. (2003, 2006). The overall better fit to the models is indicated by 

THE asterisk IN THE 'AGE' COLUMN. 



Object 


Ay [mag] 


Age [Myr] 


Mass [Mj„p] 


Te// [K] 


log g 


CAHA Tau 2 
CAHA Tau 3 
CAHA Tau 4 

CAHA Tan 5 


0.9 
4.1 
4.7 

1.0 


1* /3 

1*7 3 

1* / ^ 


8/15 
10 / 20 
15 / 50 
•■^0 / 50 


2159 / 2399 
2240 / 2637 
2337 / 2902 

2703 / 2902 


3.5 / 4.0 

3.6 / 3.9 
3.6 / 3.6 

.3.7 / 3.6 



Best fit model. 
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Av=10 N = 5251 




- J [mag] 



Fig. A1. — Color-magnitude diagram based on optical (SDSS) and NIR (Omega2000) data showing the 5251 objects fulfilling the 
selection criteria as described in the text. The NIR data is given in the CIT system, the optical filter is in the Johnson-Cousin system (see 
appendix). The solid lines with squares and diamonds are 5 Myr isochrones based on the models for low-mass stars and BDs from Baraffe 
et al. (2002) and Chabrier et al. (2000), respectively. The extinction vector is based on the interstellar extinction law of Mathis (1990) with 
Ry=3.1. The distance to Taurus is assumed to be 140 pc. 
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Fig. A2. — NIR color-color diagram for the 646 selected objects (dots) and the finally selected 16 candidates. The boxes show the 10 
already known members, while crosses indicate the positions of the final 6 low-mass candidates. The dashed lines denote the classical 
TTauri locus from Meyer ct al. (1997), the solid lines show the main sequence stars from B8 to M6 and the class III giants from GO to M7 
(Bessell & Brett 1988). The dash-dotted lines show the directions in which an M6 and an M3 dwarf would be moved due to extinction. 
The extinction law is the same as in Figure Al. The colors are given in the CIT photometric system. 
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Fig. A3. — NIR color-color diagram comparing the 6 new low-mass candidates presented in this paper (crosses), known Taurus BD from 
Guieu et al. (2006) (boxes) and embedded carbon stars from Liebert et al. (2000) (diamonds). The lines are the same as in Figure A2. The 
colors are given in the 2MASS photometric system. 
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Fig. A4. — Red/IR color-color diagrams comparing the 6 new low-mass candidates presented in this paper (crosses), known Taurus BDs 
from Guieu et al. (2006) (boxes) and very low- mass BDs in the young a Orionis cluster from Zapatero Osorio et al. (2007) (diamonds). 
The objects from Zapatero- Osorio typically have larger errors than our candidates. 



Very low-mass brown dwarfs and planetary-mass objects in Taurus 



25 




Wavelength [/^m] Wavelength [fim] 

Fig. A5. — Left: J-band spectra of CAHA Tau 1 and CAHA Tau 2 compared to those of field BDs and M-type 
(super-)giants. The reference spectra have a resolution of R ^ 2000 and were downloaded from the IRTF spectral library 
(http://irtfweb.ifa.hawaii.edu/~spex/IRTF_Spectral_Library; Gushing et al. 2005). Right: J-band spectra of CAHA Tau 1 and CAHA 
Tau 2 compared to theoretical models from Burrows et al. (2003, 2006) with R « 500 and young low-mass objects from AUers et al. (2007) 
with R « 300. Our candidate spectra were corrected with the Ay values given in Table A5 and using the extinction law of Mathis (1990) 
with Rv=3.l 
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Fig. A6. — Comparison of CAHA Tau 1 (shown 4 times in red) to low-mass field reference objects (black spectra) from Gushing et al. 
(2005). The best overall agreement is found with Kelu 1 AB a binary BD system where the combined spectrum has been established as 
reference for the spectral type L2. 
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Fig. A7. — HR-diagram comparing some of the least massive Taurus members known today. 
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Fig. A8. — Broad band SED of CAHA Tau 1 overplotted on theoretical model for a 3 Myr 10 Jupiter mass object with Te//=2040 K 
(Burrows et al. 2003, 2006). 
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Fig. A9. — NIR color-magnitude diagram (in the 2MASS system) comparing the 6 new low-mass candidates presented in this paper 
(crosses) with two sets of carbon stars. The objects from Liebert et al. (2000) lie outside the galactic plane, have assumed distances between 
~9 and ~110 kpc and are embedded in dusty shells. Nearby carbon stars of types C-J and C-N from TanaJsa et al. (2007) are shown with 
diamonds. 
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Fig. AlO. — Absolute proper motions based on a simulation using the Besangon model. Black dots are dwarfs, larger green dots are 
sub-giants, largest blue dots are giants. The candidates' relative proper motions are indicated. The grey area gives the proper motion locus 
of most Taurus members (see also Figure All). 
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Fig. All. — Proper motions of our candidates compared to proper motions of Taurus members from the literature: the compilation 
of Ducourant et al. (2005) (pale blue pluses), Hartmann et al. (1991) (green squares, objects with errors smaller than lOmas/yr in both 
directions), Frink et al. (1997) (red triangles) and Jones & Herbig (1979) (blue diamonds). Typical errors for the two latter references are 
5mas/yr. Also shown are the proper motions of stars of similar magnitudes as our two candidates CAHA Tau 1 (thick dots) and CAHA 
Tau 2 (empty circles). 
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Fig. A12. — Observed photometric values of the remaining candidates as listed in Table A5 compared to the best fit theoretical models 
shown in Table A9. The observed magnitudes were corrected for extinction using the extinction law of Mathis (1990) with Ry=3.1 and 
the estimated values for the objects' extinction given in Table A5. 



